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ABSTRACT 



Context. The process that heats the solar chromosphere is a difficult target for observational studies because the assumption of local 
thermal equilibrium (LTE) is not valid in the upper solar atmosphere, complicating the analysis of spectra. 

Aims. We investigate the linear correlation coefficient between the intensities at different wavelengths in photospheric and chromo- 
spheric spectral lines because the correlation can be determined directly for any spectra from observations or modeling. Waves that 
propagate vertically through the stratified solar atmosphere affect different wavelengths at different times when the contribution func- 
tions for each wavelength peak in different layers. This leads to a characteristic pattern of (non-)coherence of the intensity at various 
wavelengths with respect to each other that carries information on the physical processes. 

Methods. We derived the correlation matrices for several photospheric (line blends of Call H near 396 nm, 630 nm, 1082.7 nm, 
1565 nm) and chromospheric spectral lines (Call H 396.85 nm, Call IR 854 nm, Call IR 866 nm, He I 1083 nm) from observations. 
We separated locations with significant photospheric polarization signal, and thus, magnetic fields, from those without polarization 
signal. For comparison with the observations, we calculate correlation matrices for spectra from simplified LTE modeling approaches, 
1-D NLTE simulations, and a 3-D MHD simulation run. We apply the correlation method also to temperature maps at different optical 
depth layers derived from a LTE inversion of Ca II H spectra. 

Results. We find that all photospheric spectral lines show a similar pattern: a pronounced asymmetry of the correlation between line 
core and red or blue wing. The pattern cannot be reproduced with a simulation of the granulation pattern, but with waves that travel 
upwards through the formation heights of the lines. The correct asymmetry between red and blue wing only appears when a temper- 
ature enhancement occurs simultaneously with a downflow velocity in the wave simulation. All chromospheric spectral lines show a 
more complex pattern. In the case of Call H, the 1-D NLTE simulations of monochromatic waves produce a correlation matrix that 
qualitatively matches to the observations. The photospheric signature is well reproduced in the matrix derived from the 3-D MHD 
simulation. 

Conclusions. The correlation matrices of observed photospheric and chromospheric spectral lines are highly structured with char- 
acteristic and different patterns in every spectral line. The comparison with matrices derived from simulations and simple modeling 
suggests that the main driver of the detected patterns are upwards propagating waves. Application of the correlation method to 3-D 
temperature cubes seems to be a promising tool for a detailed comparison of simulation results and observations in future studies. 

Key words. Sun: chromosphere, Sun: oscillations 



1. Introduction 

The Ca II H and K spectral lines have been one of the most im- 
portant chromospheric diagnostics in solar and stellar physics. 
The lines show a generally complex behavior with some shared 
characteristics. Both H and K exhibit sudden "bright grains" 
(BGs) where wavelengths in and near the line core revert t o 
emission for a short duration (e.g. jRutten & Uitenbroeklfl99ll) . 
The B Gs are ofte n repetitive with a cadence around 180-200 sec 
(e.g.. iBeck et al.ll2008l and references therein). As explanation 
for the BGs the (upwards) propagation of acoustic waves has 
been suggested soon after their first detection, because the BGs 
can often been seen appearing in the line wing already around 
50-100 sec before the emission in the core (see Fig. [T] or iLiul 
U974tlBeck et al .1120081) . which supports the idea of a wave prop- 
agation. 

Due to their formation above the photosphere in a layer with 
mostly NLTE conditions, the theoretical treatment and the in- 
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Fig. 1. Example of observed Call H spectra. 



terpretation of the Ca line behavior in observations is difficult. 
The assumption of LTE breaks down for the formation of the 
line cores of the Call resonance lines, because of the low gas 
density, but is still valid in the wings of these lines. 

For the analysis of the spectra fortunately also diagnostic 
methods exist that are independent of the LTE or non-LTE as- 
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sumption. One of these methods is the matrix of the linear cor- 
relation coefficient of the intensity at different wavelengths. The 
correlation shows a characteristic pattern if there is a causal re- 
lationship between different atmosphere layers. When a time lag 
is introduced in the analysis of a temporal sequence of spectra, 
the pattern changes accordingly, as an intensity variation at a 
wavelength A\ at a time t\ appears at a different time h in A^. 
One of the great advantages of the correlation coefficient is that 
it can be calculated for observed or synthetic spectra from either 
numerical or analytical models in the same way, for any kind 
of spectral line or also continuum levels. Interestingly, the cor- 
relation matrices can also be derived for other physical quanti- 
ties, like temperature cuts at geometrical heights in simulations, 
which widens the field of possible applications of the analysis 
method. 

Since there is no current literature on the theoretical expecta- 
tions or t he interpretat ion of the correlation matrices (but see for 
example [Craml (Il978l) for a similar type of study in the Fourier 
domain), we want only to present in this contribution some ex- 
amples of correlation matrices for observations and four numer- 
ical experiments with a different degree of sophistication. The 
equation for calculating the correlation matrices is explained in 
Sect. [2 We show the correlation matrices for observations of 
spectral lines from near-UV (396 nm) to near-IR (1.5 /mi), cov- 
ering several photospheric and chromospheric spectral lines, in 
Sect. [3] The corresponding results for the numerical experiments 
are shown in Sect. 01 Section \5\ shows correlation matrices for 
temperature as an example of another physical quantity to which 
the method can be applied. We summarize and discuss the find- 
ings in Sect. [6l whereas Sect. [7] gives any conclusions that can 
be drawn at this time. 

2. Calculation of correlation matrices 

We use the standard definition for the linear correlation coeffi- 
cient, tab for two quantities A and B: 

r AB = ; 1 = f (A-(A)).(B-(B)), (1) 

where the integration j is to be executed over all elements of A 
and B, respectively. () denotes the average. 

For our purpose of wavelength correlations, we normally 
used a monochromatic image in one wavelength, I(Aq), as 2-D 
image A, to be correlated with B = I(A\) to obtain r^\ x . We also 
selected in some cases subfields of observations or simulations 
to investigate the influence of magnetic fields on the correla- 
tion matrices. The correlation matrices relating spectra inside the 
same spectral range are square and symmetric to the diagonal; 
i.e., r^x x = ^Mo ( see f° r example Fig. [2]). For cross-correlations 
of different spectral ranges, the resulting matrix is in general 
rectangular without a symmetry axis. Cross-correlation matri- 
ces between different spectral ranges are especially interesting 
where time-series of spectra are available, because they allow to 
study the variation of the correlation for different time lags. For 
propagating waves and a finite difference in formation height, 
the maximum correlation should be reached for a time lag corre- 
sponding to the wave travel time between the two height layers. 
Cross-correlation and correlations with a time lag will be inves- 
tigated in a subsequent publication. 

In the present case, we used both large-area scans yielding 
monochromatic images I (A, x, y) with spatial coordinates (x, y), 
and spatio-temporal time series giving I (A, x, t). Equation (Q} is 



insensitive to this difference. We note that I(A, x, y) and I(A, x, t) 
are equivalent only in a statistical sense if they both cover a suf- 
ficient large number of periods, or more general, variation times. 
In a time series, all phases of the variations are seen at a fixed 
location, whereas for a large area map the phases are sampled at 
different locations in a random state. This requires that the spa- 
tial maps should be as large as possible to avoid contamination 
of the correlation matrices by a coherent evolution in small sub- 
fields (» 5"), and that the time-series should at least cover some 
variation cycles (» 5 min). . 

3. Correlation matrices from observations 

To obtain correlation matrices for the various spectral lines, we 
chose several observations of quiet Sun (QS) areas at disc centre 
taken between 2006 and 2009. The QS data were taken with slit- 
spectrograph systems, where the 1-D slit is stepped across the 
solar image to obtain a 2-D field of view (FOV). Different spa- 
tial locations in the scanning direction are thus sampled at dif- 
ferent times. The obser vation data were a cquired with the spec- 
tropolarimeters POLIS dBeck et al.ll200l 396 nm, 630nm), TIP 
(Collados etaDl2007l 1083 nm, 1565 nm), and a spectroscopic 
setup using the main spectrograph of the German Vacuum Tower 
Telescope (VTT) for Ca II IR 866 nm. Simultaneously with the 
Ca II IR 866 nm line, also Ca II H intensity spectra were recorded 
for the investigation of cross-correlations; these spectra were, 
however, not used here. For Call IR 854 mn, we have two data 
sets. The first was recorded at the VTT with a similar setup as 
for Call IR 866 nm. The second d ata set of Ca II IR 854 nm was 
taken with the IBIS spectrometer dCaval lini 2006) of the Dunn 
Solar Telescope (Sac Peak/NSO) in spectropolarimetric mode. 
IBIS is a Fabry-Perot-Interferometer based 2-D spectrometer. 
We selected observations that covered large FOVs on or near 
disc centre, and one time-series; the Call IR 854 nm observa- 
tion with IBIS is the only exemption ( no QS, but a pore off disc 
centre). Appendix lAl shows the observed FOVs; more details on 
each observation are given in Table I A. ll 

POLIS, TIP, and IBIS measured the Stokes vector in mag- 
netic sensitive spectral lines, which additionally gives the op- 
tion to localize magnetic fields, and to avoid or include the re- 
spective positions in the analysis. To investigate the influence 
of photospheric magnetic fields on the correlation matrices, we 
created masks of magnetic field locations for each data set. 
We integrated the absolute circular polarization signal, |V|, in 
wavelength, and then set a variable threshold depending on the 
integration time and magnetic sensitivity of the observed spec- 
tral lines. Locations with polarization signal below (above) the 
threshold were assumed to be nearly field-free (magnetic) and 
are labeled "unpolarized/unpol." ("polarized/pol.") in the follow- 
ing. As third sample, we also calculated the correlation matrix 
using the full FOV, which is better suited for comparisons when 
no mask of field locations can be made. Appendix lAl shows also 
the Stokes V maps and the masks of each observation. For the 
observations of Call IR at 866 nm, no polarimetric measure- 
ments were available. We did not try to define a mask of loca- 
tions with presumably magnetic fields from the line core inten- 
sity map (see Fig. IA.3lh but only calculated the correlation matrix 
of the full FOV in that case. For the observations of Ca II IR at 
854 nm taken at the VTT in 2009, we have not yet aligned the 
data to the simultaneous polarimetric observations with TIP and 
POLIS, thus also no mask is available. 

Figures [2 [3] and |4] show the correlation matrices, either ob- 
tained by averaging the matrices of several observations in each 
wavelength range, whenever multiple suited data sets were avail- 
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Fig. 2. Same as Fig.[3]for 630nm, 1083 nm and 1565 nm. Top and bottom row of the magnification show 630.25 and 630.15 (630), 
He 1 1083 nm and Si 1 1082.7 nm (1083), and 1565.2 and 1564.8 (1565), respectively. The wavelengths for the left graphs are ,1-630 
nm, /l-1083nm and /i-1564.85nm, respectively. 



able, or derived from individual observations (Call IR lines). 
Each figure shows the correlation matrix for the full spectral 
range covered in each observation, and additionally a magnified 
view of individual spectral lines for the TIP and POLIS data sets. 
Besides for Ca II IR 866 nm and the QS data set of Ca II IR 854 
nm, each time the three samples full FOV, unpolarized and po- 
larized are shown left to right in each plot; for Call IR 866 nm 
and the QS data set of Call IR 854 nm, only the matrix of the 
full FOV is shown ( top and middle panel of Fig. [4]). We did not 
find significant changes of the correlation matrices with differ- 
ent integration times for photospheric lines, even if for example 
for 1565 nm integration times between 5 and 30 sees were used 
in the various observations (cf. Table I A. ft . The spectra a t 1.56 
micron taken with the TIP II camera (Collados et al ] l2007l) have 
a pattern of spectral fringes (see, e.g jBeck & Rez aei 2009) that 
produces high correlation on stripes parallel to the diagonal that 
are not of solar origin. We first calculated correlation matrices 
for a time- series with fixed slit separately, but they were qual- 
itatively identical to those from large-area maps and thus were 
simply included in the averaging. The He I 1083 nm line shows 
nearly no absorption in QS conditions and thus yields little struc- 
ture, but has been displayed for completeness. 



The correlation matrices for the intensity at various wave- 
lengths and spectral lines forming in photospheric atmosphere 
layers (continuum intensity, 630.15 and 630.25 nm, Si I 1082.7 
nm, 1564.8 and 1565.2 nm, the line blends in the wing of Call 
H) are very similar regardless of wavelength in near-UV, visible 
or near-IR. Correlations between the intensities at continuum 
wavelengths are always close to 1. Photospheric spectral lines 
show a reduced correlation up to anti-correlation, if the inten- 
sities at wavelengths in or close to the line core are correlated 
with the intensities at continuum wavelengths. Wavelengths in 
the red and blue wing are not equivalent: the reduction of corre- 
lation between core and red wing is always stronger than for the 
blue wing. The absolute value of the correlation is smaller for 
1565.2 nm and 630.25 nm than for 1564.8 nm or 630.15 nm. We 
ascribe this fact to the line depth of the respective spectral line. 

The correlation of the line-core intensity of a photospheric 
line with the intensity at close-by wavelengths (r^+Ai) gives 
roughly the same pattern as the correlation with the intensity 
near another line core (r^+hj) each time A2 corresponds to a 
photospheric line-core wavelength, doubling the pattern of the 
stripes of reduced correlation for each spectral line present in- 
side the wavelength range (Fig. [2]). The correlation changes only 
slightly from line to line, which is due to the fact that the re- 
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Fig. 3. Correlation matrices of the 396 nm range. Top: full ob- 
served wavelength range. Bottom graph: magnification of line 
cores of some blends (Fel at 396.455 nm, Til at 396.427 nm, 
and CrI at 396.369 nm) {top) and Call H 396.85 nm {bottom). 
The matrices of the three samples made from the FOV (full FOV, 
unpol(arized) and pol(arized)) are aligned left to right in each 
plot. A line profile is overplotted as black line. Wavelengths are 
/1-396.85 nm. White color corresponds to a correlation value of 
1. 
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Fig. 4. Top: Correlation matrix for Call IR 866 nm for the full 
FOV. Wavelengths are /i-866.215nm. Middle: correlation ma- 
trix for Call IR 854.2 nm for the full FOV (slit spectrograph 
data, disc centre). Bottom: correlation matrices for Ca II IR 854.2 
nm (IBIS data, off centre in an active region). Wavelengths are 
,1-854.215 nm. 



spective photospheric line pairs in the visible and near-IR wave- 
length ranges (1564.8/ 1565.2 nm, 630.15/630.25 nm) have sim- 
ilar formation heights (Cabrera Solanaetal .2005). If the forma- 
tion height differs strongly, like for Si I 1082.7 n and He I 1083 
nm or the line core and all the line blends of Call H, r^+Ai 
and rx x A 2 ±kA show little similarity. 

With respect to the different samples made in each FOV (full 
FOV, unpol., pol.) only small changes are seen. A clear trend 
is that for the polarized sample the width of the stripes with 
reduced correlation related to the line cores is slightly larger, 
clearly seen only for 630.25 nm and 1564.8 nm. This is due to 
the splitting of the lines inside magnetic fields, leading to multi- 
ple spectral components in the intensity profil^l The correlation 
values between line core(s) intensities and the intensities at 
continuum wavelengths are always higher in the polarized sam- 
ple than in the unpolarized sample; this is best seen for Call 
H and Si 1 1082.7 nm where the graphs of "unpol" and "pol" are 
markedly different in the value of the correlation coefficient. The 
higher correlation could be due to the field- strength dependent 
shift of the optical depth scale. The intensity will be enhanced at 
all wavelengths on locations with magnetic fields, which itself 
will not increase the correlation due to the subtraction of the av- 
erage values in Eq. (Q]). But if additionally a dependence of the 
intensity enhancement on the field strength exists, I {At, x,y) and 
I{Aj, x,y) will scale in the same way. For example, for all loca- 
tions {x, y) with stronger (weaker) than average magnetic fields, 
dh = I{A h x,y) - (I{Ai)) and dlj = I{Aj,x,y) - (I{Aj)) will be 
greater (smaller) than zero at the same time, increasing the frac- 
tion of locations with a positive correlation (dl{ • dlj > 0). 

1 Tested on synthetic spectra. 



The structures in the correlation matrices of all chromo- 
spheric Call lines are significantly different from the photo- 
spheric cases. For Call H (Fig.0, the correlation of the inten- 
sities at wavelengths in or close to the Ca line core with that at 
wavelengths separated by around 0.3 nm gives a reduced corre- 
lation over an extended wavelength range; the reduction of cor- 
relation is not restricted to only a small stripe near the line core 
like for the photospheric lines. Around the diagonal, a band of 
high correlation exists, whose width changes with wavelength, 
increasing from the Call H line core until around 0.2 nm wave- 
length separation, decreasing from 0.2 nm until 0.3 nm, and then 
increasing again. The effect is most pronounced for the polar- 
ized case, where the high correlation extends over almost the 
full wavelength range for wavelengths outside of spectral lines. 

The photospheric line blends in the Ca line wing show the 
same pattern with the red/blue asymmetry near their line core as 
the other photospheric lines, but the Call H line core shows a 
stronger structured pattern. It consists of two stripes of reduced 
correlation starting from the line core like for the photospheric 
lines, but at around 0.15 nm from the line core a small area of 
reduced correlation in the shape of a half-circle is visible that 
is missing for the other lines (white arrow in Fig. [3]). The fea- 
ture connects the two stripes of reduced correlation going to 
red and blue wing, respectively. Centred on around ±0.02 nm 
from the line core, two small squares of high correlation can be 
seen (black arrows in Fig. [3]), directly above (left) of the lower 
half-circle. The blue/red asymmetry of the photospheric lines is 
missing completely, both wings seem to be equivalent for the 
Call H line core. Call IR at 866 nm {top panel of Fig© also 
shows a more complex pattern compared to photospheric lines. 
It exhibits two patches of reduced correlation at around +0.03 
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Fig. 5. Simulation of granulation and resulting correlation matrices. 1st column: continuum intensity I c and LOS line-core velocity 
maps; tick marks are in arcsec. 2nd to 4th column: correlation matrices for 1564.8 nm, 630.15 nm, and the full Calcium line profile. 
Top row: velocity directly proportional to I c . Bottom row: same as before, but with a random velocity added. 



nm from the core, presumably caused by a strong Fe I line blend 
at this wavelength, but blue and red wing are equivalent without 
asymmetr)Q. The correlation matrix for the QS data of Ca II IR 
at 854 nm (middle panel of Fig. [4]) is fairly similar to that of the 
866 nm line, but shows only a single patch of reduced correla- 
tion. The structure seen in the off centre observations of Ca II IR 
at 854 nm is quite different. The stripes of reduced correlation 
to the red and blue of the line core do not connect directly like 
in the QS data; especially for the polarized sample a high corre- 
lation is found between the branches from core to blue and red. 
We suggest that the strong variation of the correlation matrix for 
one and the same spectral line is due to both the solar structure 
of the observed FOV (cf. Fig. IA.3b and the off centre position, 
similar to the differences in the temperature correlations later 
on in Sect. [51 In total, the photospheric lines present a simple 
structure in the correlation matrices dominated by the red/blue 
asymmetry, whereas all chromospheric lines show a more de- 
tailed fine- structure. 



4. Numerical and (semi-)analytical spectra 

For comparison to the observed correlation matrices, we used 
spectra obtained from four approaches differing in the sophisti- 
cation of the method used for the generation of the spectra. 



4.1. Granulation simulation 

The convective energy transport and the resulting granulation 
pattern dominates the spatial distribution of the continuum in- 
tensity in the solar photosphere, and also all layers up to around 
300-400 km above log rsoonm = that are relevant for the forma- 
tion of photospheric spectral lines. To estimate the contribution 
of granulation to the wavelength correlation matrices, we cre- 
ated a synthetic data set that includes only granulation, and no 
waves of any kind. We used a subsection of the continuum inten- 
sity map of a long-exposed quiet Sun observation on disc centre 
obtai ned with the TIP instrument at 1 .56 micron dBeck & Rezaeil 
2009) to define the intensity pattern in a 25 arcs 2 area (see left 



column of Fig. O. We defined the LOS velocity to be inverse 
proportional to the intensity by 



vlos (x, y) = - A/ c O, y) • 20000 ms 



-l 



(2) 



Only the wavelength range near the line core is shown in Fig. [4] 



which yields velocities in around a ± 1 kms -1 range, since 
AI c (x,y) = (I c (x,y) - (I c )) lies between around 0.95 and 1.05 
for the near-IR observations. 

The temperature stratification was defined using a modi- 
fied version of the Harvard Smithson ian Reference Atmosphere 
model (HSRA. lGingerich et al.ll 1 97 ll) . We used an optical depth 
range from log r=1.4 to log r=-6. Up to log r=-4, we used 
the values of the original HRSA model, but in the layers above 
we substituted the chromospheric temperature rise in the HSRA 
with the co rresponding values of the H olweger-Mueller model 
fflQLMUL jHolweger&Muellerll974l) . The Holweger-Mueller 
model gives to first order an atmosphere in radiative equilibrium 
without a chromospheric temperature rise. The main motivation 
for this choice was that otherwise in a LTE synthesis Ca spec- 
tra with (unobs erved) large s ingle- peaked emission in the line 
core result, but iRezaei et aD (120081) also showed that the solar 
chromosphere partially reverts to such a low-energy state in the 
absence of shock waves. The temperature stratifications of the 
two models were smoothly connected by adding the difference 
Ari ogT =_4 = 7XHSRA)i ogT =_4 - r(HOLMUL)i ogT =_4 to the val- 
ues of the HOLMUL model. The gas and electron density above 
log r=-4 were derived by an extrapolation of the exponential 
decrease of the HSRA model below log r=-4 to the rest of the 
optical depth range. The atmosphere was not put to hydrostatic 
equilibrium, as only a qualitative guess was desired. 

For the granulation simulation, the temperature on each point 
(x, y) was additionally modified by multiplying the temperature 
stratification with 1+A/ C (x, y), yielding something like ±300 K 
variation at continuum forming layers. Given temperature strat- 
ification and LOS velocity, we synthesized spectra for the 1.56 
micron range, 630 nm range, and Call H, respectively , in LTE 
with the SIR code dRuiz Cobo & del Toro Iniestalll992h . The re- 
sulting correlation matrices for 1564.8 nm, 630.15 nm, and the 
Call H line are shown in the top row of Fig. [5] from left to right. 
They show stripes of, in this case, not only reduced but anti- 
correlation starting at the line cores like in the observations, but 
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ity. Bottom: comparison of observed intensity and LOS velocity. 



do not exhibit the red/blue asymmetry. The correlation values are 
much lower than for the observations and reach down to around 
-0.9 near the very line core of any photospheric line. For Ca II H, 
the correlation basically is unity over the complete wavelength 
range outside the photospheric line blends. The asymmetry for 
photospheric lines thus does not originate from the presence of 
hot upflows and cold downflows, if a perfect linear relation be- 
tween intensity and velocity is assumed. 

As second test, we used the same setup as above for the tem- 
perature and velocity, but added a random Gaussian variation 
of 150 ms" 1 variance to the velocities derived from Eq. ([2]), and 
calculated again the correlation matrices from the corresponding 
synthetic spectra (lower row in Fig.0. The correlation matrices 
then show the a red/blue asymmetry for the photospheric lines, 
albeit the correlation values are still much lower than for the ob- 
servations. The matrix of Ca II H does not change noticeably and 
has no resemblance to the observed correlation matrix. 

Figure [6] compares the velocities of the two synthetic ap- 
proaches with the actually observed line-core velocity of the 
1564.8 nm line. The scatter of the observed velocity (blue) 
around the linear relationship (red) defined by Eq. ® is about 
three times as large as the tested 150 ms -1 variance (black). This 
suggest that the asymmetry between red and blue wing for wave- 
lengths near the line core of photospheric lines is not due to the 
granulation pattern of hot upflows and cold downflows, but the 
deviation from the purely convective source of the velocity field, 
i.e., the (magneto-)acoustic waves that usually have no strong 
signature in intensity at photospheric levels. The observed cor- 
relation matrix of the chromospheric Ca II H seems to be fully 
unrelated to a pure granulation pattern. 

Another possible source of the asymmetry in the correla- 
tion matrix of photospheric lines could be the temperature (or 
intensity) contrast inversion between granules and intergranular 
lanes above the p hotosphere ("reverse granulation", see, e.g., 
iRutten et"aT]|2004l) . To find out whether this effect influences the 
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Fig. 7. Relation between continuum intensity (I c ) and line core 
intensity (I CO re)- Top, left to right: I c , I core , polarization degree. 
Red contours trace high polarization degree, blue contours low 
I c . Middle: magnified view of the white rectangle marked above. 
Tick marks are in arcsec. Bottom: scatterplot of I c and I core . 
Red line: granulation simulation, blue dots: observation, black 
crosses: same after binning (see text for details). 



correlation matrix, we investigated the relation between contin- 
uum intensity I c and the line-core intensity I core of the 1564.8 nm 
line. Figure [7] shows the corresponding maps for the full FOV of 
the quiet Sun observations, together with a map of the polariza- 
tion degree of the spectra. The line-core intensity was defined as 
the minimum intensity value in the spectral region of the 1564.8 
nm l ine. The line-core intensity of the low-forming 1564.8 nm 
line (iCabrera Solana et al.l2005l) still reflects the granulation pat- 
tern with a small positive correlation coefficient of 0.26 with re- 
spect to 7 C , i.e., I core is not yet influenced by the reverse granula- 
tion pattern. Mainly where the polarization degree is high, I core 
and I c exhibit inverse patterns. The relation between I core and I c 
shows a steeper slope for the granulation simulation than for the 
observations (bottom panel of Fig. [7]), but this could be due to 
the stray light contribution to the observed sp ectra. The binned 
values were derived like in iBeck et all (120071) by averaging the 
pairs (I c , Icore) over bins in I c , where the plotted points corre- 
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Fig. 8. Top left: scaling law for the amplitudes of perturbations 
with optical depth. Top right: some of the temperature stratifi- 
cations for the 20 K run. Bottom: velocity perturbation for three 
different temperature perturbations. 



spond to ((F c ), (I l core )). They correspond to the center of gravity 
of the I core distribution as function of I c . From the positive corre- 
lation of I c and I core and the correlation matrix of the granulation 
simulation, we think that we can exclude the granulation pattern 
as source of the asymmetry for photospheric line, even if, e.g., 
the lines at 630 nm form in slightly higher layers closer to the 
reverse granulation pattern than the near-IR lines. 



4.2. Simplified 1-D LTE synthesis of wave propagation 

In the next attempt to reproduce the observed correlation matri- 
ces, synthetic spectra were generated by assuming an upwards 
propagating wave that creates a temperature perturbation of the 
modified HSRA model as defined above. The perturbation was 
modeled as a Gaussian with a width of <x=0.4 units in log r, 
and an amplitude A that was scaled up with decreasing optical 
depth. We scaled the amplitude with the electron pressure that is 
related to gas density using A=Aq/ ^p2 (see top panel of Fig. [8]). 
We normalized the scaling law for the amplitude to 1 at log r=0 
to be able to prescribe exact initial values at this optical depth. 
We then synthesized spectra with the SIR code for the Gaussian 
perturbation moving through all 75 points of the optical depth 
grid (log r +1.4 to -6; sampling 0.1 units of log r), only mod- 
ifying the temperature. The middle panel of Fig. [8] shows some 
of the temperature stratifications that result for an initial pertur- 
bation of 20 K at log r=0. The propagation speed of the pertur- 
bation is constant in the optical depth scale, not in geometrical 
height, and without adding a corresponding velocity perturbation 
no Doppler shifts are induced. 
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Fig. 9. Comparison of observed (left) and synthetic 1-D LTE 
spectra without velocities (middle), and with velocities included 
(right). The three runs with 10,20,30 K initial perturbation are 
shown contiguously from bottom to top in the simulated spectra. 

Figure [9] compares a sample of observed profiles at a fixed 
spatial location over a time of around 600 sees (left) with the 
synthetic spectra resulting when propagating an initial tempera- 
ture perturbation of 10,20,30 K at log r = through the optical 
depth grid (middle , from bottom to top for 10,20,30 K). We syn- 
thesized the Call H line with all of its line blends for which 
atomic parameters were at hand. The photospheric lines at 630 
nm were also synthesized to be conform with the standard wave- 
lengths ofPOLIS. 

The observed profiles have been stretched along the y-axis to 
match the evolution "speed" of the synthetic onefl The bright- 
ening s originating in the line wing and culminating in strong 
emission in the core can be clearly traced in the both synthetic 
and observed spectra. The observed bright grains increase in in- 
tensity from the 1st to the 2nd one like for the synthetic spec- 
tra, but this actually was only by chance; the observations were 
selected before the synthetic spectra were calculated. The com- 
plex behavior of the Ca line core in the observations is of course 
completely missing in the synthetic spectra due to the lack of ve- 
locities and the LTE assumption. The synthetic spectra are fully 
symmetric with respect to the Ca rest wavelength. Using the syn- 
thetic spectra for the three emulated "waves", we calculated the 
wavelength correlation matrices in the same way as for the ob- 
servations. We remark that in this case only 3x75 spectra were 
used, opposite to the observations with several ten thousands of 
spectra, but the synthetic spectra cover all phases of the wave. 
This correlation matrix, however, turned out to be fully symmet- 
ric also for the photospheric lines, i.e., the pronounced red/blue 
asymmetry (cf. Fig. [2]) was missing. 

To improve the agreement with the observations, we then 
also introduced a line-of- sight velocity perturbation in phase 
with the temperature perturbation. The velocity amplitude was 
scaled with the same relation as the temperature perturbation 
(Fig. [8l lower panel). The velocity was chosen to be positive, i.e., 
a downflow, corresponding to a temperature increase by com- 
pression of gas during the propagation of an acoustic wave. This 
choice of direction for the velocity was necessary to produce the 
correct asymmetry for the photospheric lines, whereas for the 
Ca II H line core it enhances the red emission peak, contrary to 
what happens in observations (right panel of Fig. [9]). If the ve- 
locity was applied with a time-lag, either leading or lagging the 
temperature increase, the resulting correlation matrix turned out 
to be fully symmetric again in red and blue line wing. 

Figures [TO] and [TT] show the correlation matrices for Call H 
and 630 nm, respectively, obtained from the synthetic spectra in- 
cluding the positive velocity perturbation. They have to be com- 
pared with Figs. [2] and [3] Starting with the photospheric 630 nm 
lines, the asymmetry between red and blue wing can be repro- 
duced by the assumed positive velocity. Similar to the obser- 



3 The cadence of the observations was around 20 seconds, yielding 
fewer spectra than in the wave simulation. 






Fig. 10. Wavelength correlation matrices for Call H from the synthetic 1-D LTE spectra including a velocity perturbation. Left: full 
wavelength range. Middle: magnified view of the line core. Right, same section from observations; display range was reduced to 
-0.15 to 1 here; the lower level of -0.5 applies to the left and middle plot. 



vations, the correlation coefficient is lower, when the intensity at 
the rest wavelength of the 630. 15 nm line is correlated with other 
wavelengths than for the line core of the 630.25 nm line, which 
is presumably related to the line depth. The difference of the 
absolute values of the correlation coefficient between observa- 
tions (minimum below -0.3) and the simulation (minimum +0.3) 
could be due to the granulation pattern in the observation that 
contributes anti-correlation for line-core wavelengths (Fig. [5]). 
For Ca II H, the global shape of the observed correlation matrix 
is roughly reproduced, but with several huge differences to the 
observations. The band of high correlation around the diagonal 
widens monotonically with wavelength distance to the Ca line 
core. Near the Ca line core, the half-circle/bridge between the 
stripes with reduced correlation to red and blue wing is missing; 
directly in the Ca line core, only one square of high correlation 
exists, located exactly at the rest wavelength of the Call H core. 

The fine- structure near the very line core of Call H seen in 
the correlation matrix derived from the observations thus is not 
reproduced by the simplified wave simulation. A close match, 
however, is not to be expected in this case, since the processes 
in the upper atmosphere during the passage of a shock front are 
very complex. They include large Doppler shifts of the line due 
to upwards and downwards direct motions, together with strong 
variations of the intensity near the line core, where the flow di- 
rections in front of and behind the shock front can be opposite 
dCarlsson & Stein|[l997l) . Macroscopic mass flows in the atmo- 
sphere correspond to the correlation with a shifted line profile 
that yields the same correlation matrix, but displaced in wave- 
length, if the profile shape is not changed. To test the effect of 
macroscopic flows in addition to the simulated wave, we took the 
correlation matrix for the synthetic Ca II H spectra of the wave 
simulation, shifted it along the diagonal, and added it to the orig- 
inal un- shifted one. Shifting the correlation matrix by around 5 
kms -1 and adding it to the original matrix reproduced the fea- 
tures of the observed correlation matrix near the line core, i.e, 
the double squares of high correlation and connecting bridge be- 
tween low correlation stripes, but also doubles the pattern due to 
the photospheric blends in full contradiction to the observations. 
We thus suggest that any additional macroscopic flows should 
only be present in the upper atmosphere layers not seen by the 
photospheric blends, and that the fine- structure in the correlation 
matrix near the Ca II H core could be produced by the complex 
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Fig. 11. Wavelength correlation matrix for 630 nm from the syn- 
thetic 1-D LTE spectra including a velocity perturbation. 

flow pattern near a shock front that passes through the atmo- 
sphere. 

The two contributions from the granulation simulation and 
the simplified wave propagation together seem to be able to re- 
produce the patterns seen in the correlation matrix of observed 
photospheric spectra and the global structure of the Ca II H cor- 
relation matrix, albeit not the fine structure near the line core of 
the chromo spheric line. 

4.3. NLTE 1-D simulations 

As third numerical experiment, we used spectra generated 
from 1-D NLTE simulation s don e with the WAVE code 
(iRammacher & Ulmschneide rl |2Q03 | ). T he simulations were 
done similar to ICarlsson & Steinl (|1997|) . with a photospheric 
piston that excites acoustic waves of various types. The driver 
of the piston was varied between stochastic excitation and 
monochromatic waves. Here we use three runs corresponding 
to stochastic waves with a mechanical energy flux of 5-10 7 erg 
cm" 2 s _1 and 2 - 10 8 erg cm~ 2 s _1 , and monochromatic waves of 45 
sec period with also 2-10 8 erg cm~ 2 s _1 as energy flux. 

The correct calculation of spectra for the chromospheric 
emission lines requires the use of partial redistribution 
(PRD). However, a line treatment implementing PRD like in 
Ulmschneider (1994) leads to excessive computation times that 
cannot be tole rated in time-dependent wave c alculations. We 
thus followed iHuenerth & Ulmschneide j (1 1995b by employing 
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Fig. 12. Left to right: correlation matrices of the Call H line core from 1-D NLTE calculations for stochastic waves (1,2), monochro- 
matic waves (3), corresponding section of matrix from observations. Wavelength range (-0.03 nm to 0.03 nm) and display range 
(0.15 to 1) of the observations' matrix are slightly different. Wavelengths are /i-396.849nm. 



a so-called "pseudo-PRD", which computes the line assuming 
complete frequency redistribution (CRD) but then artificially re- 
moves the damping wings from the lines by multiplying the 
damping parameter in the Voigt function by a factor of 1/100. 
The spectra were calculated with this method from the full sim- 
ulation runs only for the very core of Ca II H in a range of ±0.03 
nm around the rest wavelength, but this covers the fine-structure 
as seen in the observations. Spectral samples of 1500 such pro- 
files were used to calculate the correlation matrices for each pis- 
ton model. 

Figure [12] shows the resulting correlation matrices for the 
runs with stochastic excitation (1 and 2), and the monochromatic 
case (3). The patterns in the correlation matrices for stochastic 
excitation do not match at all to the observations {last column 
in Fig. [12]). The matrix for the monochromatic waves, however, 
fits partly to that of the observed spectra. The central band of 
high correlation is smaller than for the observations and there 
exists a cross-like structure of higher correlation values at the 
rest wavelength, but the high correlation values at ±0.025 nm 
appear in both observations and simulations. The bridge/half- 
circle that connects the two bands with reduced correlation start- 
ing from the core is present in the simulations as well. Some of 
the additional fine- structure in the simulations' correlation ma- 
trix is missing in the observations, but this actually could also 
be due to either limitations on the observational side (tempo- 
ral/spatial/spectral resolution) or the realism of the simulations 
with a strictly monochromatic permanent driver. 

To investigate the correlation matrix for the 1-D NLTE 
simulations also in the line wings, new spectra were synthe- 
sized from the simulation res ults using the PRD code as de- 
scribed in lUlmschneiderl (1 19941) . This allows to extend the wave- 
length range, but unfortunately is quite demanding on computing 
power. So far, only a set of 100 spectra with 3 seconds temporal 
sampling are available for the case of the monochromatic 45- 
sec piston (Fig.[13j top). The spectra were calculated on a non- 
equidistant wavelength grid with a spectral sampling of around 
4 pm in the line wing and 1 .5 pm near the line core to reduce the 
computational effort. We interpolated and re- sampled the spectra 
to an equidistant sampling of 2 pm for a better match to the ob- 
servations (1.9 p m sampling at a spectral resolution of 220.000, 
iBeck et aDl2005l) . and then determined the correlation matrix as 
before. The middle and bottom panel in Fig. [13] show the re- 
sulting correlation matrix for the full wavelength range and a 
magnification of the line core, respectively. Similar to the simpli- 
fied 1-D LTE modeling, the correlation band around the diagonal 
widens monotonically for increasing distance from the line core, 
contrary to the observations (cf. Fig. [3]). The correlation from 
core to red and blue wing (ro,^. and r^o) now shows positive 
correlation for \Xi\ > 0.1 nm, contrary to both the observations 
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Fig. 13. Top to bottom: 1-D NLTE spectra for monochromatic 
waves, correlation matrix of Call H, magnification of line core. 
Wavelengths are /i-396.849nm. Display range is ±1. 



and the simplified LTE modeling. This could be due to the short 
time span sampled and the strict period of the driver; the correla- 
tion of a spectrally un- smoothed version of the spectra without 
re- sampling shows several repeated stripes of alternating high 
and low correlation there. The line core region {right) matches 
better to the observations (cp. to the rightmost panel in Fig.fTOb. 
The two squares of high correlation at ±0.02 nm are present, as 
well as the bridge between the stripes with reduced correlation 
to red and blue wing, albeit the correlation only drops to zero 
and does not turn to anti-correlation (r+0.2,+0.2). 

In total, the 1-D LTE or NLTE calculations are able to re- 
produce the general shape of the correlation matrix of Ca II H 
in the line wing, whereas the 1-D NLTE calculation yields a 
first order match to the very line core. In both cases, the lower 
part of the atmosphere is treated similarly as being permeated 
by propagating acoustic waves leading to similar spectral pat- 
terns (cp. the spectra in Figs. [9] and [13]). The line core region 
is only reproduced in the NLTE calculations due to the complex 
dynamics in the shock fronts that form in the upper atmosphere 
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Fig. 14. Correlation matrices from 3-D MHD simulations for 
630.15 nm (top left), 630.25 nm (top right), 1564.8 nm (bottom 
left), and 1565.2 nm (bottom right). 



(e.g.. iRammach er & Ulm schneiderl [l992l) . The remaining mis- 
match between observed and synthetic correlation matrices is 
presumably due to the granulation pattern that is absent in the 
synthetic spectra of the wave simulations. 

4.4. 3-D MHD simulation 

As last numerical model, we obtained spectra in the 1565 nm 
and 630 nm range fr om a 3-D MHD sim ulation run done with 
the C0 5 BOLD code dFrevtag et aTll2002h. The s i mulations ar e 
described in more detail in S chaff enberge r et all (120051 120061) ; 
the simulation box had an extension of 4.8 Mm x 4.8 Mm in 
the horizontal domain. The spectra were calculated for a sin- 
gle snapshot of a simulation run that included magnetic fields; 
the dynamics in the upper atmosphere is, however, dominated 
by dynamic events and sho ck fronts like in the field-free case 
(cf. lWedemeyer et al]|2004l) . The radiative transfer is treaded by 
mean opacities in the simulation, which prevents to obtain NLTE 
spectra of, e.g., the chromospheric Call lines directly from the 
simulation results without additional calculations; we thus only 
use the photospheric spectra of Fel 630.15 nm, 630.25, 1564.8 
nm and 1565.2 nm that can be synthesized at once. Figure [14] 
shows the correlation matrices for these four spectral lines ob- 
tained from the 3-D MHD simulation (cp. to Fig. [2]). The struc- 
ture of the correlation matrix from the observed spectra is repro- 
duced accurately: the lines with larger line depth (630.15/1564.8 
nm) show a stronger reduction of correlation, their stripes with 
reduced correlation are broader, and all lines in the matrix from 
the simulation show the correct red/blue asymmetry. The sole 
difference to the observations is the magnitude of the correla- 
tion coefficients: the simulation generally shows stronger anti- 
correlation (display range from -0.7 instead of -0.3 like for the 
observations), which presumably is due to the higher spatial res- 
olution of the simulation (40 km grid size). The correlation ma- 
trix from the 3-D MHD simulation matches much better to the 
observations than the simplified 1-D LTE modeling of wave 
propagation, again presumably due to the absence of the gran- 
ulation pattern in the latter. 



5. Correlation matrices for temperature 

To highlight the possible usefulness of the correlation matrix 
for an analysis of also other physical quantities than spec- 



tra, we also calculated correlation matrices for temperature as 
function of optical depth. The temperature values were derived 
by an inversion of C all H spectra of an 1-hour time-series 
(cf. iBeck etaD 120081 observation nr. 21 in the present paper, 
cf. Appendix lA.lb and three scans of an active region off disc 
centre (taken on 08/12/2007, 9 ~ 50°; not listed in the ap- 
pendix ). The inversion is based on the SIR code in synthe sis 
mode dRuiz Cobo & del Toro Iniestalll992l:lRuiz Coboll 19981) : it 
will be discussed in more detail in a forthcoming paper. It mainly 
differs from the standard version of SIR in the usage of a pre- 
calculated fixed intensity-temperature response function, and the 
initial use of an archive of pre-calculated Ca spectra. The usage 
of the archive was found to be necessary to achieve a reasonable 
match of the Ca line core region; the standard iterative approach 
of SIR failed there because of the complex temperature stratifi- 
cations required. The inversion with the LTE assumption basi- 
cally maps the intensity values at a given wavelength in the Ca 
profile to some range in optical de pth, like given by the intensity 
cont ribution function (cf. Fig. 5 of lRezaei et al . 2008; B eck et all 
2009, A&A, in press). 

Figure [15] shows the resulting temperature maps at six levels 
(logr=0,-l,...,-5) at left. The temperature maps of the time-series 
are to first order identical to the intensities in the spe ctral win- 
dows OW, MW2, MW3, IW1 H 2 V and Ca line core in lBeck et al.1 
(2008, Fig. 4). Correlating the maps of T(Ti, x, t) and T(tj, x, t) 
using Eq. (Q]) yielded the correlation matrices shown at the right 
in Fig. [15] For the time-series, the shape of the temperature cor- 
relation does not match well to any pattern in the Ca II H matrix, 
but the matrix of Ca II IR at 854 nm (Fig. [4]) shows some similar- 
ity. The structure seen from AA = -0.05 nm to nm in the matrix 
from the spectra matches to that from log r= -3 to -6 in tempera- 
ture: an iterative sequence of increased and reduced width of the 
high correlation band around the diagonal with lower correla- 
tion at, e.g., rAi=-o.o5,o and ri ogT =_3 _5, respectively. The matrices 
from the temperature maps of the active region are similar to 
each other, but differ from the one for quiet Sun, which will be 
due to the different heliocentric angle and the solar surface struc- 
ture like for the two observations of Ca II IR 854 nm. They only 
keep a few common properties, like, e.g., the anti-correlation at 
around logr=-2.5 to all regions above. A conversion from log r 
to the corresponding wavelength (or vice versa) will be needed 
for a direct comparison with the correlation matrices from spec- 
tra, or a conversion from T(r) to geometrical height for a com- 
parison with simulations, but the inversion code first has to be 
improved for a quantitative analysis. 

6. Summary and discussion 

We have calculated the linear correlation coefficient between in- 
tensities at different wavelengths for various photospheric and 
chromospheric spectral lines from near-UV to near-IR wave- 
lengths. The correlation coefficients yield a matrix that contains 
information on the causal relationship between different wave- 
lengths, and thus, height layers in the solar atmosphere. For all 
photospheric spectral lines, a pronounced asymmetry between 
the red and blue wing is found. The intensity at the line core 
has a weaker correlation with the intensities in the red wing than 
with those in the blue one. Correlations between the intensities 
at continuum wavelengths are always high. All chromospheric 
spectral lines considered (Call H,CaII IR at 854 nm and 866 
nm) show much more structure in the correlation matrix near 
their respective line cores than the photospheric lines. 

A simple granulation simulation that assumes hot upflows 
and cold downflows produces a strong reduction of the corre- 
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Fig. 15. Left top: temperature maps at various optical depths from a LTE inversion of Call H spectra of a time-series. Right top: 
correlation matrix of temperature as function of optical depth. Bottom three rows: same for three scans of an active off disc centre. 



lation between line-core intensities of photospheric lines and 
the intensity at continuum wavelengths. The resulting correla- 
tion matrix, however, is fully symmetric in red and blue wing. 
Synthetic LTE spectra corresponding to the upwards passage 
of a temperature perturbation are able to reproduce the addi- 
tional signature seen in the photospheric lines, the pronounced 
asymmetry between red and blue wing, if a positive velocity 
(downflow) is introduced together with the temperature enhance- 
ment. Spectra of Call H synthesized from 1-D NLTE simula- 
tions, where propagating waves are generated by a photospheric 
piston, yield a correlation matrix that qualitatively matches to 
the pattern in the correlation matrix of observed spectra near the 
very line core of Ca II H, but only for the case of a monochro- 
matic driver with 45 sec period. Stochastic excitation of the pis- 
ton produces a correlation matrix that is inconsistent with the ob- 
servations. A 3-D MHD simulation run, in which the dynamics 
in the upper atmosphere is dominated by shock waves, faithfully 
reproduces the asymmetry of red and blue wing seen for all pho- 
tospheric lines, together with the relative amount of correlation 
that seems to depend on the line depth of the respective lines. 
The 3-D simulations contain contributions from both granula- 
tion and waves. Temperature maps at different optical depth lay- 
ers that were derived from an LTE inversion of observed spectra 
yield a correlation matrix that has at least a structural resem- 
blance to the matrices obtained from observed spectra of Call 



IR 854 nm. A clear difference between quiet Sun and active re- 
gions is seen in the temperature correlation. 



Influence of granulation pattern The granulation pattern leaves 
its imprint on the observed correlation matrices in all photo- 
spheric spectral lines, which partially hides the signature of the 
dynamic processes relevant for the chromosphere. The observed 
asymmetry between the intensities of red and blue wing, how- 
ever, could not be related to the granulation pattern (Sect. 14. lb . 
To better isolate the dynamic processes, it seems advisable to 
partly remove the contribution of to the granulation pattern from 
the spectra prior to the calculation of wavelength correlation 
matrices in future studies. For large-area maps, this could be 
achieved by a local (or global) adjustment of continuum intensity 
values. One could force all spectra in a small-scale (1-3" radius) 
surrounding of each pixel to an identical intensity in continuum 
wavelengths, corresponding to removing the granulation pattern 
locally, and then calculate several "local" correlation matrices to 
be averaged later on, or force the continuum intensity to unity 
across the full FOV. Time-series will be even better suited for 
the removal of the granulation pattern, because the granulation 
can be more accurately filtered in the temporal than in the spatial 
domain: the typical time- scale of granules (5 min) differs from 
the fast evolving and propagating waves (around 1 min typical 
life time from first visibility in continuum layers to shock sig- 



12 



C.Beck & W.Rammacher: Correlation matrices of spectral lines 



nature in the chromosphere, e.g jBeck et al.ll2008|) . whereas the 
spatial scales of granules and chromospheric brightenings can 
be of comparable size. The signature of the dynamic processes 
will then become more prominent in granulation- filtered data. 
Preliminary tests with, e.g., the method of local intensity bal- 
ancing significantly reduced the areas of high correlation around 
the diagonal. 

7. Conclusions 

The matrices of the linear correlation coefficient of photospheric 
and chromospheric spectral lines carry information on the phys- 
ical processes in the solar atmosphere. They are influenced by 
both the convective granulation pattern and more transient events 
like propagating waves that are presumably related to the chro- 
mospheric heating process. We have shown that the most pro- 
nounced feature in the correlation matrices of photospheric lines, 
an asymmetry between the intensities in the red and blue wing, 
can be reproduced to first order by assuming propagating (acous- 
tic) waves that produce an intensity enhancement together with 
a red- shift of the spectral lines. Such propagating waves also 
qualitatively reproduce the correlation matrix of chromospheric 
spectral lines. 

The calculation of correlation matrices can be applied to 
other physical quantities as well, like for example temperature 
maps at different height levels in the solar atmosphere. The 
method seems suited to compare observations of chromospheric 
spectral lines, where NLTE conditions apply, with theoretical 
or numerical modeling of wave propagation in the solar atmo- 
sphere. The analysis of observed spectra can be refined by fil- 
tering out the signature of the granulation pattern that differs 
mainly in its temporal behavior from the more transient chro- 
mospheric heating process. Other physical quantities like 3-D 
temperature cubes can be analyzed with the correlation method 
as well, which seems to be an useful tool for a comparison of 
simulations and observations. 
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Fig. A.l. Examples of large-area scans used. Left to right, con- 
tinuum intensity, absolute integrated Stokes V signal, mask of 
magnetic fields. Top to bottom: 1565-6 (see. Table IA. ll for the 
details, the number -X denotes the observation); 1083-7; 630-10; 
396-15. For Call H, the integrated line core emission (H-index) 
is shown instead of V. 




Fig. A.2. Overview over the time- series used, same layout as 
Fig. ED Top to bottom: 396-21, 630-22, 1083-23. 



Appendix A: Used observations 

A.1. TIP/POLIS 

For the derivation of the correlation matrices, several data sets 
taken with the spectropolarimeters TIP and POLIS were used. 
Table IA. 1 1 lists the data sets with a consecutive number, the op- 
eration number of the data set on the day, the date, the total in- 
tegration time per slit position, and the approximate extent of 
the scanned area. For each of the polarimetric data sets (TIP at 
1083 nm or 1565 nm, POLIS at 630 nm), we show a contin- 
uum intensity map at left, the absolute wavelength integrated 
circular polarization in the middle, and the mask of field-free 
and magnetic locations derived from the polarization signal at 
right. For the TIP data in the 1083 nm range, the polarization 
signal of the photospheric Si I line at 1082.7 nm was used to 
define the masks. Some of the TIP observations consist of re- 
peated scans over the same area, we have not cut them apart for 



Table A.l. Detailed list of the observations used. The operation 
number (2nd column) is the number of the observation on that 
day. 



Nr. 


Op. 


wavelength 


ddmmyy 


integ. 


size x/y 




nr 


nm 




time 


in arc sec 


1 


3 


1565 


10/08/06 


30 sec 


35/70 


2 


3 


1565 


12/08/06 


10 


120/70 


3 


1 


1565 


14/08/06 


10 


200/70 


4 


1 


1565 


21/08/06 


6 


75/70 


5 


1 


1565 


07/09/07 


5 


80/70 


6 


7 


1565 


07/09/07 


5 


80/70 


7 


5 


1083 


24/07/06 


6.5 


75/35 


8 


6 


1083 


24/07/06 


6.5 


75/35 


9 


7 


1083 


24/07/06 


6.5 


75/35 


10 


6 


630 


24/07/06 


6.5 


75/87 


11 


7 


630 


24/07/06 


6.5 


75/87 


12 


8 


630 


24/07/06 


6.5 


75/87 


13 


2 


396 


18/07/06 


6.5 


40/63 


14 


3 


396 


18/07/06 


6.5 


40/63 


15 


6 


396 


24/07/06 


6.5 


75/63 


16 


7 


396 


24/07/06 


6.5 


75/63 


17 


9 


396 


08/12/07 


3.3 


100/63 


18 




866 


20/06/07 


0.25 


120/177 


19 





854 


23/09/08 


100 ms 1 


39/72 


20 




854 


29/08/09 


25 sec 


60/84 


21 


1 


396 


24/07/06 


3.3 


52763 


22 


1 


630 


24/07/06 


3.3 


52 2 /63 


23 


1 


1083 


24/07/06 


3 


52 2 /35 



IBIS data, exposure time per wavelength; 2 : duration in min 



the display but show them contiguously. For the Call H data 
from POLIS, the intensity integrated over 0.1 nm around the Ca 
line core (H-index) is shown in the middle. The mask in that 
case has been created from the simultaneous 630 nm data of the 
second POLIS channel. The maps for the time-series used here 
(Fig. IA.2b are organized in the same way as for the large-area 
scans. Spatial sampling in scanning direction was 075, along the 
slit the sampling was 07 15 (630 nm), a/29 (396 nm), 0717 (1083 
and 1565 nm). Additional information on the respective data 
sets can be found in their online overview archives (accessible 
from http://www.kis.uni-freiburg.de — > Observatories — > Data 
archives). 

A.2. Ca II IR 

The spectroscopic data of Call IR 866 nm was taken at the VTT 
using a PCO camera at the main spectrograph. Simultaneously 
with the Ca II IR 866 nm line also Ca II H was observed. A study 
of centre-to-limb variation and some time series are available 
from the observation campaign and can be used in further in- 
vestigations. As no polarimetric data is available in this case, 
Fig. lA.3l only shows a pseudo-continuum map of the intensity in 
the line wing, and the wavelength integrated line core intensity 
comparable to the H-index for Call H. 

For Call IR 854 nm, we used two different data sets. The 
first map (Ca II IR 854-20) was obtained again at the VTT using 
a PCO camera at the main spectrograph. The camera was run in 
parallel to TIP in He I 1083 nm and POLIS in its standard con- 
figuration; the data were taken in quiet Sun at disc centre. We 
have not yet aligned the data of all the instruments, so we have 
not been able to create a mask of locations with significant polar- 
ization signal and used the full FOV instead. The second data set 
for Call IR 854 nm was taken in September 2008 with the IBIS 





Fig. A.3. Top: overview of the Call IR map 866-18. Left: wing 
intensity, right: core intensity. Middle: overview of the Call IR 
map 854-20. Left: wing intensity, right: core intensity. Bottom: 
overview of the Call IR map 854-19. Left to right: wing inten- 
sity, core intensity, mask. 



spectrometer at the Dunn Solar Telescope in Sac Peak/NM. The 
instrument was used in the spectropolarimetric mode. The mask 
was in this case derived from the absolute integrated Stokes V 
signal of Call IR 854 nm for simplicity; a second data set in 
630.25 nm taken two minutes later is also available. The 854 
nm line was scanned in 30 steps with a spectral sampling of 4.3 
pm; spatial sampling was 07 17. The data set was on an active re- 
gion containing a pore; it is not fully compatible to the quiet Sun 
observations on disc centre that were used for all other spectral 
lines. 



